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Summary 

I-  »■>  L' 

Linear  autonomous  functional  differential  equations  of  neutral  type 
include  Volterra  integral  and  integrodif ferential  equations  as  special 
cases.  The  paper  considers  numerical  approximation  of  solutions  to 
these  equations  by  first  converting  the  initial  value  problem  to  an 
abstract  Cauchy  problem  in  a  product  space  *  weighted  L^-space)  and 
then  using  abstract  approximation  results  for  C  -semigroups  combined  with 
Galerkin  type  ideas.  In  order  to  obtain  concrete  schemes  subspaces  of 
Legendre  and  Laguerre  polynomials  are  used.  The  convergence  properties 
of  the  algorithms  are  demonstrated  by  several  examples. 


Running  head:  Volterra  type  equations 


Subject  classification:  34G10,  34K99 ,  45D05 ,  45J05.45L10 


1.  Introduction 


The  purpose  of  this  paper  is  to  develop  an  approximation  scherr. 

O 

based  on  [^-approximation  by  orthogonal  polynomials  for  the 
following  problem: 

p  t  to 

y(t)  =  <>  +  l  A.  f  x(T-h.)dT  +  J  J  A(o )x(x+o )dodx 

j=0  J  6  J  0 

t 

+  J  f(i  )dx ,  t  >  0,  (1.1) 

0 

P  o 

x(t)  =  y(t)  +  \  B.x(t-h.)  +  J  B(x)x(t+x)dx  a.e.  on  t 

j=l  J  J 

x  ( t  )  -  (t)  a.e.  on  t  <  0, 


where  0  =  h  <  h  <  ...  <  h  =  h,  <j>°  e  ]Rn  and  4)1  resp.  f  is 

n  u  i  p 

an  IR  -valued  function  on  (-“,0]  resp.  (0,=°).  Furthermore,  A., 

J 

B.  are  n*n-matrices  and  A(*),  B( • )  are  n*n-matrix  valued 

J 

functions  on  (-®,0  1.  It  will  be  convenient  to  define 

P  0 

D(x  )  =  x(t)  -  l  B.x(t-h.)  -  /  B(x  )x(t  +  x )dx , 

C  j=l  J  J 

p  0 

L(x^)  =  j  A.x(t-h.)  +  J  A( x )x(t  +  x  )dx  , 

1  j=0  J  J  -- 

where  as  usual  for  a  function  x:  ]R  -*•  IRn  the  functions 

xfc  :  ( -  °° ,  0  ]  -*■  IRn,  t  >  0,  are  defined  by  x^Cx)  =  x(t  +  x),  x  _<  0 . 

The  state  of  problem  (1.1)  at  time  t  naturally  is  the  pair 

(y (t ) ,x  ) .  Correspondingly  we  choose  as  a  state  space 

Z  =  IRn  x  L^(-°°,0;  IRn)  with  norm 
g 

I  4>  I  2  '  U°|2  +  J  U1  ( T  )  I  2g( T  )dx  ,  <t>  =  U0,^1)  e  Z, 

—  ao 

where  the  weighting  function  g  is  of  the  form 


,  .  Bx 

g( x )  =  e  ,  x  <  0  , 


(1.2) 


with  >  0  , 


•  j  is  the  Euclidean  norm  on  1R  . 


Before  we  discuss  the  solution  semigroup  of  problem  (1.1)  in  the 
next  section  we  indicate  some  special  cases  covered  by  (1.1). 

It  is  clear  that  y(t;<t>)  is  always  continuous  on  t  >_  0 , 

whereas  x(t;<j>)  need  not  to  be  continuous  on  t  _>  0 .  This  motivates 

to  introduce  the  pair  (y(t),xf.)  as  the  state  of  (1.1)  at  time  t 

and  not  the  pair  (x(t),x  )  (see  (  3  1).  If  for  a  solution  of  (1.1) 

x(t;<j>)  is  continuous  on  ]R  then  y(t)  =  D(x  )  for  all  t  _>  0, 

0  1  t 

4>  =  D  ( ),  y(t)  is  locally  absolutely  continuous  a.e.  on  t  _>  C 

and  x(t)  is  a  solution  of 


-^r  D(x  )  =  L(x  )  +  f^t)  a.e.  on  t  >  0, 

CA  V  \y  L/ 

x(t)  =  <t'1(t)  for  all  t  <  0. 


(1.3) 


Equation  (1.3)  is  a  functional-differential  equation  of  neutral 

type.  Of  course,  if  E.  =  0,  j  =  l,...,p,  and  B  =  0  then  (1.3) 

J 

is  an  equation  of  retarded  type.  Solutions  of  (1.1)  for  general 
4>  C:  Z  could  be  considered  as  generalized  solutions  of  (1.3).  In 
this  case  y(t)  =  D  ( xfc  )  only  a.e.  on  t  _>  0 .  Further  important 
types  of  equations  covered  by  (1.1)  are  Volterra  int  egro-di  f  ferer.t  i  a 
and  integral  equations, 

t 

x ( t )  =  J  A  (t-T )x(t  )dx  ,  t  >  0 

0  (1.-4) 

x  ( 0  )  =  4>°  , 


(flj  - 

II 

CQ 

O 

o 

CXJ 

III 

0, 

A(o )  =  A ^ (-o  )  , 

♦l  ^ 

0) 

and 

t 

x  ( t ) 

=  f(t)  + 

J 

B  ( t -t  ) X ( T  )d T  , 

C+ 

1  v 

0, 

0 

<AJ  = 

°;  Bj  = 

111 

< 

•\ 

o 

0, 

B( o )  =  B1 (-o  )  , 

♦l  = 

0, 

f(t)  = 

whe  re 

f  is 

locally 

absolutely  continuous 

on 

1  V 

O 

(1.5) 


f  (  T  )d  7  ), 


2.  The  solution  semigroup 

We  impose  the  following  condition: 

A,  B  c-  L^/g(  — ,0;  IRnXn).  (2.1) 

Since  the  weighting  function  g  defined  in  (1.2)  trivially 
satisfies  hypotheses  (HI)  and  (H2)  of  [  9  ]  we  get  immediately 
from  !  9  ;  Thm  2.11. 

Proposition  2.1.  Assume  f  =  0.  Then  the  family  T(t),  t  _>  0 , 
of  operators  defined  by 

T  ( t  )  t>  -  ( y  ( t  ;  <i )  ,  xt  ( <J> ) )  ,  t  _>  0,  <j>  e  Z  , 

where  y(t;$),  x(t;4>)  is  the  solution  of  (1.1)  corresponding  to 
<t> ,  is  a  CQ-semigroup  on  Z. 

Let  A  be  the  infinitesimal  generator  of  T(*).  Then 

Proposition  2.2.  A  is  given  by 

0  1  1 

dom  A  =  { ( <{>  ,  $  )  e  Z  |  <j>  is  locally  absolutely  continuous 
on  (--,01,  i1  e  L^(--,0;lRn)  and  DU1)  =  <j>° }  , 

A(<t>°U1)  =  (L( -ti1  )  U1  )  ,  ( <}>0 , 1  )  e  dom  A. 

Proof :  Let  <t>  e  D(A)  and  choose  a  e  TR  sufficiently  large.  Then 

oo 

”1  —  X  t 

<t>  =  (AI-A)  tp  -  J  e  T(t)ipdt  for  a  p  e  Z  which  is  equivalent  to 

0 

4>  =  /  e  y  (t  )dt , 

0 

00  oo 

41(t)  =  |  e  Atx(t  +  T  )dt  =  e*T  J  e 


(  2 . 2  ) 


■At 


x  ( t  ;  )  d  t  ,  t  <  0  .  ( 2  .  3  ) 


Equation  (2.3)  shows  that  c  is  locally  absolutely  continuous 
on  (-<*  ,0  ]  and 


"1  1  *1  o 

♦  =  -  <P  €  L„(— ,0  ;  IRn )  .  (2.1 j  ) 

C 

1  11 

From  A <f»  -  [  A<p  ]  -  o  and  (2.^4)  we  see 

1  •  1 
[A<H  =  <P  . 

Taking  Laplace-transforms  in  the  second  equation  of  (1.1)  and 
observing  (2.2),  (2.3)  we  get 

<fr°  =  DU1). 

Differentiating  the  first  equation  in  (1.1)  and  then  taking 
Laplace-transforms  we  obtain  analogously 

A$°  -  =  LU1). 

This  and  A  <}> 0  -  [  A4>  )  ^  shows 

0  1 

[A*]  =  LU  ). 


Thus  we  have  shown  that  the  operator  given  in  the  proposition 
is  an  extension  of  the  infinitesimal  generator  of  T(-).  Call 
this  extension  for  the  moment  C  and  choose  <t>  e  dom  C.  We  put, 
for  A  sufficiently  large,  ip  -  (  A I  —  C  )  <J>  and  4>  =  (  AI  —  A)-1^  e  dor.  A, 

i.e.  =  (  A I  —  A )  <(>  .  Then 

A 

p°  =  A  D  (  4> 1  )  -  L  U 1  )  =  ADUb  -  L(*h, 


These  two  equations  imply 


<t>a(T)  -  <ty  (  t  )  "  eAl(iJi1(0)  -  4>|(0)), 


and 

[  AD  ( e  A  *  I  )  -  L(eX‘l)  ](4»1(0)  -  *J(0))  =  0. 


The  estimate 


|I  -  D  (  e  A "  I  )  +  y  L(e  A  *  I )  | 


p  -  An . 

£  l  !  B  .  1  e 

~  j=l 


J 


J  *  (rbs)1/2|Bi,2  *  i  I 


Jl/g 


j=0 


'V* 


-Ah . 
J 


1,1  ,1/2  |  .  | 

+  A  2  A  +  6  A  T  2 

1/g 

A  •  A  * 

shows  that  [AD(e  I)  -  L(e  I)]  exists  for  A  sufficiently 

1  1 

large,  which  implies  <j>  (0)  =  <J>*(0)  and  therefore  also  $  =  <t>  e  dorr.  A 


It  will  be  necessary  to  consider  problem  (1.1)  in  state  spaces 
with  weighting  functions  different  from  g  as  defined  in  (1.2).  Let 
0  <  y  <  2 , 


g(x)-eYT,  t^O,  (2.5) 

and  put  Z  =  !Rn  *  L?(-»,0;IRn)  with  the  usual  norm.  Then 

obviously  A,  B  e  L?  (-00  ;0  ;  ]Rn  n)  and  therefore  Propositions  2.1 

1 '  & 

and  2.2  are  also  true  for  problem  (1.1)  considered  in  the  state 
space  Z.  Let  T(>)  and  A  denote  the  solution  semigroup  of  (1.1) 
in  Z  and  its  infinitesimal  generator,  respectively. 

Lemma  2.3-  a)  Z  is  dense  in  Z  and  the  embedding  Z  c  Z  is 
cont i nuous . 


t)  rot’  any  v  >  Lt(-°°,0;  P  )  and  any  a  >  0  we  have 

i  i  i  av  f  l9.(-°°,0  ;  Pn). 
b 

c)  For  any  v  e  Lh(-“,0;IRn)  such  that  also  v  c  L' ( -°°  ,0  ;  IF.n  ) 

6  S 

any  a  >  0  we  have 


lin  e"  i 

i  a  ,  . 

T  |  v  (  X  ) 

T  -►  —  OO 

et  s  D, 
k 

:  =  dom  A 

"  k  k 

d)  The  sets  :=  dom  A'  c  dom  A  ,  k  =  1,2,...,  are  dense  in  Z 

and  invariant  with  respect  to  T(*). 

Proo  f :  a)  Density  of  Z  follows  from  density  of  L' ( -<*  , 0  ;  pn  )  i 

2  n  ^ 

L  (_oo,0;lR  ).  The  latter  property  is  obvious,  because  all  continu 

functions  with  compact  support  are  contained  in  Z.  Since  y  <  g, 

we  obviously  have  \  <*.  i  n  <  i$ld  for  all  $  e  Z. 

b)  This  follows  from 


, ■  ■ «  ,2 
T  I  V  ;  . 


r  ( P~ Y  )  x  .  [2a  y  T  |  ,  ,  |  2 
Je  .  t  !  e  I  v  (  t  )  !  d  x 


2  a  .  2  a  |  1 2  T  2  /  „ 

- )  1  v  ,  ,  v  e  L-(-°°,0; 

(6-y)e  If'  g 


pn) . 


c)  Using  v ( t  )  =  v ( 0  )  +  /  v(x)dT  we  obtain 

-t  °-i 

e 2  |  t  |  a  |  v ( t ) |  <  e^  j  t  |  a | v ( 0 ) | 


+  e  i  t  |  { j  e  do) 

6  T 

^  e2  j  x  |  a  |  v  ( 0  )  |  +  f- 

which  implies  the  result, 
ci  )  Clearly,  D^_  is  dense  in  Z  and  invariant  with  respect  to 
T(*)  (cf.  for  instance  1141). 


2T|  I  a . 9  - y  c  .1/2,9  Y  o  |  •  ,  .  ,  2  .1/2 
i  t  j  ( J  e  do)  ( J  e  1  v ( o  )  ,  do) 


2  (  P  "  Y  ) T 


X  !  !  v  ;  ?  ,  T  2.  o 
if' 


From  a)  it  follows  that  D  is  dense  in  Z.  Since  T(t)  Z  = 
ir  is  clear  that  D  is  also  invariant  with  respect  to  T(* 
The  inclusion  c  is  oivious  from  a)  and  Proposition  2 

For  the  rionhomopeneous  equation  we  have 

Proposition  Let  f  e  L^q  (0 ,  °° ;  ]Rn )  and  let  x(t),  y(t 

the  solution  of  (1.1).  Then 

t 

(y(t),x  )  =  T  ( t  )  <j>  +  /  T(t-x  )  (  f  (  t  )  ,0  )di  ,  t  >  0. 

0 

For  equations  with  finite  delays  (2.6)  was  proved  in  I  3  ] 
(see  also  [  15  ] ,  Section  2.3)-  The  proof  for  the  infinite  de 
case  is  quite  analogous  and  is  left  to  the  reader. 


}.  Legendre  and  Laguerre  polynomials 

In  this  section  we  state  repectively  prove  convergence  results 
2 

concerning  L  -approximation  by  Legendre  and  Laguerre  polynomia-5. 


3.1.  Legendre  polynomials 

For  the  convenience  of  the  reader  we  first  collect  some  well- 
known  facts  on  Legendre  polynomials  (see  for  instance  11]).  The 
n-th  Legendre  polynomial  Pn(t),  n  =  0,1,...,  is  of  degree  n  s 
satisfies 

[  ( 1 -t 2  )  P  '  ]  +  n(n+l)P  =  0,  P  (1)  =  1. 
dt  n  n  1  n 


For  all  n  -  0,1,... 


Pn<1)  = 
lpn(t)l 


Pn(1)  = 


1,  P„(-l)  =  (-l)n, 

1  1  on  [-1,1], 

n(n  +  l )  A  ,  1 v  r-1  ninill 
2  *  nv  ;  k  '  2 


( ?  • : : 


The  sequence  P  (t),  n  =  0,1,..,  forms  a  complete  orthogonal 
2  n 

set  in  L  ( -1  , 1  ;  IR  )  , 


1 

J  F  (t)P  (t  )dt 
i  n  m 


'  tibi  for  n  ~~ 

0  for  n  i  m. 


(3.2) 


The  derivative  of  P  (t)  is  a  polynomial  of  degree 
a  combination  of  P^ (t ) , . . . ,Pn_^ (t ) , 


P2k(t  ) 


kI 1  (^j-o)P2i  +  1(t), 
j  -  0 


p 

2k*  1 


(t  ) 


k 

=  l  (4j  +  l)P  .(t), 
j=0  J 


n  -  1  and  t  hu 


( 


1 

J 


k 


0,1,2,.. 


Using  the  orthogonality  relations  (3-2)  we  get  for 
IF ( -  1  , 1  ;  F  )  the  expansion 


f  =  i  f . p . ,  f .  =  lii! 
j=0  J  J  J  2 


/  f(t)P.(t)dt. 

2  -1  J 


he  the  orthogonal  projection  of  f  onto  span ( P  , . . . ,Pr > 
If  ;  -  1  , 1  ;  F. )  ,  i  .  e  . 


I  r,p, 


■  r  h:e  following  convergence  results.  When  there  is  no 
r.f  Ion  L'  denotes  di  f  ferentiation. 


-  0,1,...  there  exists  a  constart 


,,k  ’ 

W r  ’ '  -1,1;  IP. ) . 

I,,,...  there  exists  a  constant  c  =  c(k)  such 


<  - _  1  f*  i 

: '  _  L-  _  1  /  2  1  1  1  L-  p 

N  K  J  w K  ’  d 


WK  ’M-l  ,1  ;  F) 


any  k  -  1,2,...  there  exists  a  constant  c  =  c(k)  such  that 


r(-D  -  f  en;  <  ^rT72lrlwk;2 


for  all  r  £  w  ’  (-1  ,1  ;  TO. 


Proof.  a)  and  b)  are  special  versions  of  results  obtained  in  [  t  ] . 
c)  is  given  in  [  7  ]  ■ 


0,1,2,...  , 


is  of  degree  n 


3-2.  Laguerre  polynomials 

The  Laguerre  polynomial  L^,  n  - 

sat  is  fies 

tL"  +  (l-t)L'  +  nL  =0,  L  (0)  =  1.  (  3 . L 

n  n  n  ’  n 


The  sequence  of  Laguerre  polynomials  Ln,L  , ...  is  complete  and 

2  u  l 

orthogonal  in  L  (0,“;  1R)  with  weight  w(t)  -  e  ,  t  >  0, 


oo 


J 

0 


e  fcL  (t )L  (t )dt 
m  n 


1  for  m  =  n, 

0  for  m  /  n. 


(  3  • : 


For  the  derivatives  we  have  in  analogy  to  (3.3)  the  formula 

d  n_1 

£  L  =  -  I  L.,  n  =  0,1,...,  (3.6 

d.  n  J=0  j 

-1 

(as  usual  £  =  0  ) . 

j=0 

In  order  to  derive  convergence  results  analogous  to  those  of 
Theorem  3.1  we  need  some  preparation.  On  the  linear  subspaces 

8^  =  {v  e  C2k  ^(0,°°;  IR  |  v^2k  1  ^  locally  absolutely  conti 

on  (0,-),  tmv(j)  e  L2  ( 0  ,  °° :  IP. )  ,  m  =  0,1,...,  j  =  0,.. 

and  lim  e  o/,2tmv^^(t)  =  0,  m  =  0,1,...,  j  =  0,..., 
t  -+-00 

k  =  1,2,...,  we  define  the  operator  B  by 

(Bv)(t)  =  tv(t)  +  (l-t)v(t),  v  e  B  . 

Lemma  3-2.  a)  L  e  8,  and  L  =  BL  for  n  =  0,1,2,..., 
- - -  n  k  n  n  n  *  *  *  * 

k  -  1,2,...  .  Moreover 


Bl  v  c  8^,  j  =  0,...,k-l,  for  all  v  e  8^. 


L'  )  Let,  v  c  C  (0 ;  IR ) .  Then  v  e  8  implies  v  e  8  ,  k  =  1, 

2  K 

c)  B  is  symmetric  in  LT(0,®;3R). 

Proo f .  The  first  part  of  a)  is  trivial  because  of  (3.4).  Let 
v  e  6 ,  .  Then  v  is  a  linear  combination  of  terms  of  the  fo 
tuv*'v  ,  u  =  0,1,2,...,  v  =  0,1,.  . . ,2j  ,  and  therefore  tm(3“v) 
trn(BJv)"  are  linear  combinations  of  terms  of  the  form  twv^v 
m  =  0,1,...,  v  =  0,l,...,2j+2  £  2k.  Then  the  result  is  obviou 
b  )We  only  have  to  prove  tmv  e  L2(0,°°;1R),  m  =  0,1,...,  and 

_ ^  /  2  T1  ^ 

e  tv(t)  ->-0  as  t  -+  ® ,  m  =  0,1,...  .  Since  trivially 

tmv(0)  e  L2(0,®;1R)  and  e  t//^tIT'v(0)  ->-0  as  t  -*•  °°  we  only 
w  t- 

m  ^  * 

to  investigate  t  J  v(t  )dx .  The  result  then  follows  from  the 


est imates 


J  e  tt2m  ||  v(T)dT|2dt  1  J  e  tt2m+1  J  |v(-r)|2dTdt 


=  J  e  1 | v ( t  )  |  p(i  )dx  <  ®, 

0 

where  p(t)  is  a  polynomial  of  degree  2m+l,  and 

:  -t/2.m  \  •,  ..  I  -t/2.m+l/2t*| 

I e  t  J  v(T  )dx |  <  e  t  I v |  2  • 

0  L 

w 

c)  Density  of  8^  is  clear  by  a).  Let  u~,v  e  8^.  Then 


e  ( 1-t  )u(t  )v(t  ) 


-t/2  ...  ~t/2  ...  -t/2.  ...  -t/2  ...  . 

=  e  u(t)e  v(t)  -  e  tu(t)e  v(t)->0  as  t 


and  also 


te  Lu(t)v(t)  =  e  t  ^  2t u ( t  )e  u/^v(t)  -*•  0, 


-t/2 


te  u(t)v(t)-*-0  as  t  -*■  ®. 


Therefore  by  integration  by  parts 


<Bu,v>  3  -  |  e  t[tu(t)  +  (1-t  )u(t )  ]v(t  )dt 


w 


=  -  u(0)v(0)  -  j  e  t((l-t)v(t)  +  tv(t)]u(t)dt 
0 


-  J  e  t[(l-t)v(t)  -  (2-t )v (t ) ] u(t  )dt 
0 


=  J  e  t((l-t)v(t)  +  tv(t))u(t)dt  =  <u,Bv> 


0 


w 


We  now  are  in  a  position  to  prove  convergence  results  for 
Laguerre  polynomials  similar  to  those  of  Theorem  3.1  for  Legendre 
polynomials  following  the  approach  taken  in  t  4  ].  For 


u  =  l  uRLk  e  L^(0  uR  =  J  e_t  u  ( t  )L  ( t  )dt  ,  k  =  0,1,2.... 

k  =  0  „  0 

let  u1^  be  the  image  of  u  under  the  orthogonal  projection 


PK  f  (0  ;  IR  )  -*■  span(LQ  , .  .  .L  ) ,  i  ,e. 


N 


N  dN  r  T 

u  =  P  u  =  l  u  L  . 

k-0 


Theorem  3 . 3 •  Lot  u  e  B  .  Then 


N,  1 

u  -  u  |  p  <_  - 

L  (N+l) 
w 


£  1 Bku |  2  ,  N  =  0,1,... 
Lw 


Proof .  Using  Lemma  3.2  we  get 


u .  -  <u , L . >  „ 

J  J  l2 

w 


1  1 
—  <u,BL.>  _  =  -  —  <Bu,L.>  ~ 
J  J  l2  J  J  l2 

w  w 


_/  1 \ k  Dk  T 

...  -  ( —  "r )  <  B  u  ,  L  .  >  p. 

J  J  yt 

w 


This  implies 


}« 


L  j=N  +  l 
w  ° 

1 


.  1  „  .2k 

J=M+1  J 


<B  u  ,L. >  ? 
J  L 

w 


1  ,ak  1 2 

2  k  '  ®  u  '  ?  * 
(N+l)  L 

w 


N  =  0,1,... 


It  will  be  convenient  to  use  the  notations  D  =  -rr  and 

N  N  dt 

=  P  D  -  DP  . 


Lemma  3-^-  a)  For  any  u  e  span ( LQ , . . . ,L^ ) 


|Du i  ?  <  N | u 1  ? ,  N  =  0 ,1 ,2 , . . . 


b)  For  any  u  such  that  u  e  8, 

k 

|KNu|  p  <  ~A,y  |Bkul  p  ,  N  =  1,2,...  . 
if  nk  1  ^  i/ 

w  w 

N 

Proof .  a)  We  have  u  =  \  u.L.  and  using  (3.6) 

j=0  J  J 

N  N  j-1  N-l  N 

Du  =  l  U.L.  =  -  lu.  y  L.  =  -  l  (  y  u.)L 

j=l  J  J  j--r  J  i  =0  1  i  =  0  j=i  +  l  J 

Therefore  by  Cauchy’s  inequality 

-  N-l  N  _  N-l  N  N 

lDu|  2  =  l  <  l  u,)2!  I  (  l  1)(  I  |u,  f) 
L  i  =  0  j  =i  + 1  J  i =0  j  = i  + 1  j  = i  + 1  J 


N-l  N 

l  (N-i)  l 


j  =  i  + 1 


l  Iu,|2  l  (N-i) 


j=0  j-1 

oo 

b)  Let  u  =  y  z.L..  Then 

*-  l  l 


,  /  v  * \  v  i  1 2  N ( N  + 1  )  i  1 2  m2  i  1 2 

1  (  l  _  J  )  L  I  u  j  |  <  — p I  u  |  p  <  N  |  u  |  . 


N  V- 

P  Du  =  y  z .L. . 

j=0  J  J 


We  next  compute  the  Fourier-coefficients  of  J  u(o)do. 

0 

u  e  6  we  have  for  some  constant  c  >  0 

K 


j  u ( x  )  |  <  ce  ,t>0. 


This  implies 


/  e  T | L - ( x  )  |  /  [ u ( a  )  | d a  <  2c  /  e  T | L . ( t  )  |  ( e  -1  )d t 
0  J  0  0  J 
GO 

<  2c  j  e  T//^|L.(T)dx  <  <*>. 

0  J 

Therefore  we  can  use  Fubini ' s . theorem  in  the  following  comoutatio 

J 

For  the  moment  we  put  w.  -  I  z.L.  and  get 

J  i  =  0  11 

OO  X  CX)  oo 

j  e  tL.(t)  J  u(o)dodi  =  J  u(o)  J  e  TL.(x)dTdo 
0  J  0  0  o  J 

=  J  e  °u(c)(  l  L^(o))do  =  <u,  £  L-^>  ? 

0  i =0  J  i=0  J  L 

w 

=  <w.,  I  lF '>  ~  -  I  e  °w.(o)(  y  lF  Fdo 

J  1=0  J  L2  0  J  1=0  J 

w 


J  e  tL.(t)  /  w . (o )dodi . 


From  (3-6)  we  obtain 


/  L^(o)do  =  ( t )  -  L^+^(t) 


and  therefore 


T  .  J  T  J 

J  w.(o)do  =  £  z  .  J  L.(o)do  =  £  z . ( L. ( t )  -  L.  ( i ) ) 
0  J  i  -0  1  0  1  i=0  1  1  1 


---in-'  this  above  we  get 


00  _T  T  j 

J  e  tL.(t)  J  u(o)dodx  =  )  z .  <L .  , 


L.-L,  .  .  > 


i  =  0  1  J”  1  1+1  L2 


Zj  '  Zj-1  f°r  J'  =  ^2,..., 
zn  for  j  -  0. 


From  u(t)  =  u(0)  +  J  u(o)do  we  obtain 

0 

00 

u  =(u(0)  +  z  )L  +  £  (z.-z.  .)L. 

00  1-1  J  J“1  J 


PM  N 

u  z  p  Li  =  ( u  ( 0  )  +  2  0  )  L„  +  l  (z.-z.  )L 

j  -  ^  J  J  j 

Using,  (3-6)  this  implies 


N-l  N 


Dp  “  =  I  <=  -2  >  l  L  -  l  (  l  <Z  -Z  >) 

1=1  J1  J  A  -  n  1  1  -n  i -i  I  J1  J 


j=l  J"1  J  i  =0 
N-l 

=  ,  L(zi'ZM)Li‘ 


i=0  j  =i  +  l 


This  and  (3-8)  give 


K  u  =  l  zjLj  "  l  (zi"ZN)Li 
j=0  J  J  j -0  J  N  J 

N-l  N 

-  z..l  +  y  zml .  =  z.,  y  l . . 

N  N  j--o  N  J  N  j=o  J 

Using  Theorem  3.3  for  N-l  and  u  we  get 

I  KNu  I  2  --  (N  +  l  )  |  z  !  2  <  ( N+  1  )  l  |  z  ,  |  2  =  (N+l  )  |  u 


_  N  +  l  lnk  -  ,2  ...  2  ink  -  ,2 

-  2 k  u I  2  -  '  2k-l  I B  u I  2 

N  L  N  1  L 

w  w 


Theorem  3-5-  a)  Let  u  be  such  that  ;  e  8  .  Then  the 

„  ,  _  k 


following  estimate  is  true: 


lr-,L  N.  |  c 

!D  <u_u  >!  2  -  0^171/2  '  N  = 

Lw  N 

for  l  =  1,2,...,  where  c  =  c(;Bku|  , . . . , | Bku ( * ^  j  „). 

Lw  L2 

.  w  w 

b)  Let  u  be  such  that  u  e  B  .  Then  for  any  T  >  0 

K 


L  (0  ,T;  1R)  N 


k-1/2  ’ 


where  c  =  c  (T,  j  B  u  \  ?  ) . 

'  L2 
w 

c)  Let  u  be  such  that  u  e  6  .  Then 

K 

[u(0)  -  (PNu)(0)|  <  i  j  Bku |  2,  N  =  1,2,...  . 

Lw 

Proof .  a)  An  easy  induction  shows 

D*  -  D*PN  --  D*  -  PV  +  Y  D^V-1^',  £  =  1,2,...  . 

j=0 

Using  t.  hii  s  formula.  Theorem  3-3  for  u  ^ 1  ^  and  Lemma  J> .  H  (note, 
that  K  u  e  span ( L^ , . . . , L^ ) )  we  get 

|D*(u-uN)|  <  — - - r-  |BkuU)| 

Lw<1M>  L 


*  T  nJ  ~t£t2  lBku('-j)|  2 

J  =0  NK  Uci  L2 

w 

<  -J1 _  (  lo*  <j>, 

-  k-£  +  l  /2  1  B  u  2  * 

N  j  =  1  L 


j  =0  N 


k-1/2 


b)  Since  lul  2  <  eT  2 | u |  _ ,  the  result  is  an  immediate 

L  (0  ,T;  IR)  L2 

w 

consequence  of  part  a)  for  £  M  and  Sobolev's  embedding  theorem 


$ 

':l 


c)  Let  u  =  £  u.L.  and  u  =  £  z.L..  Then  for  J  -  1,?,... 


u.  =  <u,L.>  _  =  -  —  <u,BL.>  „ 
j  ’  J  l2  j  J  L2 

w  w 


Us  ing 


we  get 


1  r  d  -t  d 

j  q  dt  6  dt 


6  at  lj  = 


u  .  =  J  e  t ( L  .  - L . _  )udt  = 
J  n  J  <-  1 


00 

L.)u(t)dt  =  —  j  te  (-^7  L. 

J  J  A  dt  j 


)u(t )dt 


z  .  -  z  .  . 

J  J-l 


Then 


u  ( 0  )  -  (  P  u  )  ( 0  )  =  u  ( 0  )  - 


=  u(0)  - 


=  u  ( 0  )  - 


ZN ' 


I  ui  =  u(0)  -  u  -  l  (z  .-z  .  1  ) 
j=0  J  U  j  =  l  J  J  ‘ 


U0  +  z0  ZN 


j  e  tu(t)dt  +  /  e  fcu(t)dt  -  z. 
0  0 


From  u  e  8k  we  get 


I  k* 

z .  -  -  — r  <  B  u  ,  L  .  >  0  ,  J  =  1,2,...  . 
J  ,-k  Jr2 

J  Lw 


There  fore 


!  u  ( 0  )  -  (PNu)(0)|2  =  |  zN  |  2  <  l  |z.|2 

N  j=N  J 


V  _JL_  |  Rk  •  r 
l  pu-  I  <B  u 


r2  1  j  Rk  •  I  2 

>  <  ~^rr:  B  u  , 


j  zN  j" 


j  r  2  N2k 


=  1,2, 


-  18  - 


^ .  Formulation  of  the  approximation  scheme 

Following  the  general  outline  given  in  [12]  or  [111  v;e  form 
the  approximation  scheme  for  problem  (1.1)  using  Legendre  and 
Laguerre  polynomials.  Throughout  this  section  we  replace  the 
weighting  function  g  as  defined  in  (1.2)  by 


g(  i  ) 


f  1  for  -h  <  t  <  0 } 


.  eB(T+h)  for  t  <  -h 


(- 


where  as  before  6  >  C.  This  weighting  function  gives  the  same 

space  with  an  equivalent  norm  as  the  one  defined  in  (1.2).  Let 

r.  =  h.-h.  i  =  1 , .  .  .  ,p,  and  define  the  functions  £, .  , 
iii-l  i 

i  =  1  , . . .  ,p  +  l ,  by 


for 


i 


;.(t)  -  F~(hi-1+hi+2T ) *  "hi  -  T 

i 

=  1 , . . . ,p  and  by 
C  a(t)  =  -  6 ( t  +  h ) ,  T  <  -h. 


_hi-l* 


(*» 


(~ 


Trivial  computations  show  that  the  mappings  0^  defined 
0 ^  x  =  x  o  ,  i  =  1  , .  .  .  ,p  + 1 ,  are  metric  isomorphisms 

L2(-l,l;F>n)  -  L2(-h.,-h.  ;Fn),  i  =  l,...,p,  and  L2(C 

2  n  li-l  w 

Lq(-°°,-h;  1R  ),  respectively. 


by 


»;  3F.n 


We  put 


f  P.U.(t))I  for  -h.  <  t  <  -h.  , 

li  l  —  l-l 


eij(T)  =  < 


(4 


0  elsewhere , 


l 


0  , . . .  ,N  ,  and 


elsewhere 


j  =  0,...,N.  Furthermore,  let 


e00  =  =  (0,eij),  i  =  l,...,p+l,  j  =  0,...,N, 

E'  =  ( 6  1 0  ’  ‘  *  ‘  ,ep+!  5  E  =  (eoo  ,610  ’  '  ‘  '  ,ep+l  ,N  ) 


Y.  =  span(e.0 , . . . ,e-N)  ,  i  =  l,...,p+l, 

x  x 

Y  Y1  *  "  p  +  l  * 

Z"  z  ^  *  yN  :  sPan(e00’ei0’*--’ep  +  l,NK 

N  N 

Let  p  :  Z  ■*  Z  be  the  orthogonal  projection.  The  coordinate 
N  N  0  1 

vector  a  ( p  $ )  of  4  -  ( <t>  ,  $  )  e  Z  ic  given  by 


N  N  N  - 1  0  1 

a  (p  <t>)  -  (Q  )  col(h  5<e01  ><T  >  p 

L 


'p+1  ,  N 


whe  re 


„n  ;n  ;n  , .  /T  n  n  l  T  s  , ,,  , . 

Q  =  <E  ,L  >z  =  diag(I,r1q  ,...,r  q  In(N  +  1))»  ^-b) 

where  Pn(p+i)  denotes  the  n(N+l)  *  n(N+l )-identity  matrix  and 

*n-‘  r-<<eij’eik>I2)j,k=o,...,N  13  given  by 

1  Lr. 

N  ^  .  1  1  x  „  T 

q  -  diag(  1  ,y  .  .  .  »2N4l  )  ®  I- 


For  later  use  we  note  that  for  <t>  ,  ^  e  Z  with  coordinate  vector 
N  N 

a  (it),  «  (  iJj  )  ,  respectively,  we  have 


Following  the  general  scheme  as  outlined  in  (121  or  111] 


.N 


we  introduce  the  approximating  de lta- impulses  6^  by 


„N  „N  N  U  .  n  ,T 

oi  -  E  Yj,  Yi  =  (Q  )  E  (-h.-0)  , 


i  =  0  , .  . .  ,p.  Easy  computations  using  (  4  .  4)  ,  (4.5)  and  ( 4 .  6  )  she 
t  hat 


M  -1 

Y  /  -  — —  col  (C  ,...,0,1,3,...  ,2N+1  ,0  , .  .  .  ,0  )  ®  I,  i  =  0 


i  + 1 


where  the  nonzero  entries  occur  at  positions  i  to  i  +  N,  and 


Y  ’  =  8  col ( 0  ,  .  . .  ,0 , 1  , . .  .  ,  1 )  ®  I . 

r 


N 


For  later  use  we  compute  the  norm  of  6^  considered  as  an  operate 

]Rn  -*■  Ll  .  Using  (4.7)  and  the  explicit  representation  for  y^  give 

g  '  n  1 

above  we  get  for  x  <=  IF. 


6 Vx  | 


XN  JI  T.  N  sT^N  N 

<6.x,c.x>  2  -  X  ( Yi  )  Q  Yix 


N 


r1 —  l  C2k+ 1 )  i x | 2 
i  +  1  k=0 


(N+l)2  ,  | 2 

-  x 


ri  + 1 


for  i  =  0,...,p-l,  i.e. 


,N  /  1  •  r\  a 

6. |  =  (- )  (N+l),  i  =  0 , . . . ,p- 1 

i  +  1 


(4.8) 


Analogous  computations  give 


> . N i  , N+l, 1/2 

'V  =  (~} 


6 


(4.9) 


21 


A in  i  i  (see  [roof  of  Lemma  5.1)  we  obtain 

<6.'x,«f'L>  =  x1  p*  (-h.  -0  )  ,  i  =  0  , .  .  .  ,p  ,  {< 

1  L" 

€ 

Pi  1  *1 

for  x  r  Hi  and  $  e  y  .  in  analogy  to  the  definition  giver. 

in  ’  12  1  for  the  retarded  finite  delay  case  we  define  the 

•  „  •  „  ,  N 

approximating  operators  cy 

(  A'*  i>  )  ^  =  A „  ( t  ^  f  1 .  i  1  ( - h .  )  +  f  B  (  t  )  i  ^  (  t  )  d  i  ) 


•  _*  i 


A  .  <t 1  (  -h  .  )  +  J  A  (  t  )  (J>  ( t  ) d  ■ 


(  A  .  1 1 


A rvC  +  2  (A.+AnB.  )<t1(-h.  )  +  /[A(i)  +  AnE(-r)]^1(T)dT, 

U  •  _  yy  1  U  1  1  U 


(A^)1  =  JV  +  4l(^°-i1(0)+  l  B.(f1(-h.  )+  J  B(t  )<f1(T  )dx  ) 

G  T  U  •  _  *  1  1 

1  “°°  (A.  12) 


6VU  (-h.  )  -  $  (-fr -0  )  )  , 


'or  <t>  =  (t0,*1)  e  ZU. 


The  approximation  z  (t)  to  (y(t),x  )  is  given  by 

zN(t)  =  e A  t  p N  <J>  +  /  eA  (t_s  ^  (  f  (s  )  ,0  )ds  ,  t  >  0,  (^.13) 


7,f,(t)  =  ANzK(t)  +  (f(t),0),  t  >  0, 

N  N 

z  ( 0  )  =  p  if . 


For  the  implementation  of  the  scheme  we  have  to  compute  matrix 

N  N 

representations  [A  ]  for  the-  operators  A  .  As  in  [12]  we  get 


■’/.v v  >v-v  v-i 


„*v  *.  /.  ,%  .  .  .  •.  . 

-  -  ^  „  *-  J  *  -  *  h  1  »  '  .  M  •* 

 «  b.  k.  K  n  4  4  a  H  M  M  - 


ifnAiLlijUl  ikJLkJl  »i  »■  .jl  ti 


where 


HN  =  <en,an£nv 


Using  the  definition  of  the  basis  elements  and  ( H .  1 1  )  ,  (^.12)  we 
get 


Ao  =  (V60>’ 


ANe,  ,  =  ( (-1 )^  +  1(A.  +  A  B.  )  +  A.  .  +  A  B .  .  ,  ^  e.  . 

i  o  i  it  n  1 1 5  t 


!J  '  '  ' '  '■  1  o  i '  "ij  ■  “0~ij  5  dT  cij 

+  l-  --  ,  .  v j  + 1 . N  rN 


.  N , 


+  «0((-D,J  Bi  +  Bij)  +  (“DJ  <5^  - 
i  -  ^  •  jP  )  j  -  05,„.jN5 


A  ep+1.J  '  (VlJ  +  A0Bp+l,j’  di  ep+l,j  +  60Bp+l ,j 


-  <5 


j  =  0 , . . . ,N ,  where 


-h.  _ 

J1  A(x)e.  (x)dx  for  i  =  l,...,p, 
-h  ^ 


/  A(t  )eD  +  i  -(i)dx  for  i  =  p  +  1, 

Oo  *  5 


and 


-h  .  _ 

\  J1  B(x)e.  .  (x)dx  for  i  =  l,...,p. 


f  i  net;  -r—  e .  .  is  a  polynomial  of  decree  k-1  on  (-h.  ,-h.  „  ) 

QTlK  "  i  i  - 1 

for  i  =  l,...,p  and  on  (-°°,-h)  for  i  =  p  +  1  ,  respectively, 
wo  ret  immediately  from  the  orthogonality  relations  (3.2)  ar.d  ( 


<eij’  dT  eik%2  =  0  for  J  i  k. 


t  case  j  <  k  we  get 


:eir  It  eik>T2  =  l1  ei ,(T)4.k(T)dT  =  (i-(-DJ+Ji 

L  -  h . 

e  1 

-  <eik>lr  ei  j  >.2  - 

e  < 


i  =  and  similarly 


<eP+l,j’  dx  eP+l,k>L2  "  I* 


Using  (4.10)  and  the  definition  of  the  basis  elements  we  obtain 


<e00  ,A  e00>Z  =  A0 ’ 


<^00,A  ^  i  j >  Z  =  ('1)J  {Ai  +  A0Bi}  +  Aij  +  Vij’ 


i  =  1  ,  .  .  .  ,p  ,  j  =  0  , .  .  .  ,N  , 


<e00*A  ep  +  l,j>z  =  Ap+l,j  +  A0Bp  +  l,j’  J’  = 

N  - 

<o^j,A  ®qq>2  ~  ^  >  j  “  0,...,N, 

N-  (i+  if 

<e  .  , A  e  >  -  <e  .  —  e  >  —  I  +  f-  iih  +  B 

lj  '  lk  Z  lj  *  dx  Ik  l2  1  11  *lk’ 

6 

=  0 , .  .  .  ,N  , 


<e.  ;  , ANe  .  ,  =  (-1  )kB,  +  B.  ,  , 


j,k  =  0,...,N,  i  =  All  other  inner  products  are  zero. 

From  this  and  (^.15),  (A.l6)  it  follows  that  is  given  by 


wr.e  re 


4  -  (Ai*fl0Bi*fli0*A0Bi0’'<Ai*A0Bi>+Ail*fl0Bil>--' ’("1)N(Ai' 


+  A .  ,  +  A  B .  , 
lN  0  l! 


i  -  lj...,p5 


ap+l  =  ( AP  +  1 ,0  +  A0Bp+l  ,0  ’ '  '  ■  ,AP+1 ,N+A0Bp+1 ,N) 


y 


X 


course,  the  subspaces  Y.  need  not  to  have  the  same 

1  N 

dimension,  1  =  l,...,p+l.  We  could  also  take  Y^  =  span(e.  , . . .  ,e . 
i  -  The  resulting  modifications  are  obvious. 

ft  T 

Of  course,  it  is  also  possible  to  choose  g(0  =e  ,  g  >  0, 
and  to  defined  ;(  t  )  =  -8t,  t  <  0.  Then  0\  =  x  o  t;  defines  an 
isomorphism  L2(0,“;lRn)  -+  L2  ( - » ,0  ;  IRn )  .  The  basis  elements  of  Z" 

w  g 

are  defined  by 


e  .  (  t 


-  Lj  U  (t)  )I,  t  <  0  ,  j  =  0,  .  .  .  ,N, 


e00  r  ej  -  j  =  0 , .  .  .  , N . 

Of  course,  Y1J  =  span  ( eQ  , .  .  .  ,e )  and  ZN  =  ]Rn  *  yN.  In  this  case 
we  have 

=  diag(I  ,qK) , 

where  ah  -  -  diag( I , .  .  .  ,1 )  e  ]Rn ( N+ 1 ) x n ( N+ 1  ^ .  We  onlv  have  to 
introduce  6^  =  (en  , .  .  .  ,eN )  Yq  with  yJJ  =  6  col  ( I  , .  .  .  ,1  ) .  Analogou 


- —  0  o 

to  (4.9)  we  eet 


.  / N+l si/ 2 

t0  1  "  (  8  ^ 


The  operators  A  are  defined  by 


(AN0)°  =  Aq  4>°  +  \  (Ai  +  A0B.)$1(-h.)+  /  [  A(  t  )  + AqB(  t  )  ]$1(i)di 


(4.18) 


(AN«)1  =  4> 1  +  6  ( 4>°- 4> 1  ( 0  )  +  l  B.4>1(-h.)+  Jb  (  t  )  m 1  (  t  )  d  i ) 
u  i  =  l  1  1  -» 

for  $  -  <  ZN.  Again  ( AN  ]  is  given  by  (4.14).  By  analogou 

computations  as  in  the  previous  case  we  obtain 


.  tree!  01  convergence  for  the  approximation  scheme 

The  following  assumption  on  the  matrices  , . . .  ,  B  will  t< 

A)  There  exist  numbers  A.  >  o,  j  -  1  , .  .  .  ,p ,  such  that  the 

J 

np*np-matrix 


(Br..B  )  - 


is  negative  definit. 


2(1+  y  A .  ) 


diag( A  I  , .  .  .  ,  ^ 


Hypothesis  (A)  is  certainly  satisfied  if  |B.|,  j 
sufficiently  small.  Moreover,  by  a  transformation 


y  ( t ) 


'x( t  )  ,  a  >  o , 


we  car.  always  transform  system  (1.1)  into  a  system  such  that  (A 

is  satisfied.  Let  A,  B  satisfy  (2.1)  and  let  tf1  e  L2(-“,C;I- 

6  t  £ 

where  g(i)  =  e  with  some  B  e  JR.  System  (1.1)  is  equivaler. 


lx  =  o  +j  Lx  dc  +  /  f ( s )ds ,  t  >  0  a. e . , 
Z  C  S  0 

x(t)  =  0  ( t  )  ,  t  <  0  a.e. 


Then  y(t)  satisfies 


^oy  -  e  at(t^  +  /  e  S^L  yodx  +  e  at  j  f(s)ds,  *. 

0  a  S  0 


y  (t  )  =  e  4>  (t  )  ,  t  <  0  a.e.  , 


where 


p  -ah.  0 

D  y  =  y(t)  -  l  B.e  Jy(t-h.)-  f  B(t  )e°Ty(t +  t  ; It  , 

L  j=l  J  J 


9] 


I.  v  =  )  A-C  Jy(t-h.)  +  f  A(i)e  y(t  +  i)dT. 

J  ;  0  J  J  -  oo 

t  he  equation  for  y  we  get 

t  t  „  t  s  ,  V 

-  r,  ,  r  -as,  0  r  r  'O  S‘T  .  ,, 

a  i  D  v  ds  =  j  ae  ds<$>  +  J  a  j  e  L  y  did: 

.?  a"S  0  oo  a  T 

0  -at  0  \  -a  ( t-S  )  T  ,  „  ^  \  ■ 

«  -  e  ^  ~  J  e  Laysds  +  I  ' 

0  b  0 

t  t 

-  r-  at  J  f  (  t  ) d t  +  /  e  a  S  f  ( s  )  d  s  . 

0  G 

Are  v  satisfies  the  equation 

t  t 

L  v  =  $ v'  +  \  ( L  v  -  aD  y  )ds  +  (  e  a"f(s)ds,  t 

<  ■  t  .  a  ‘  s  a J  s  A 

0  0 

y  ,  *.  -  e  u  1  ;  ‘  ( t  )  ,  t  <  0  a  .  e . 


a  t  .  a  i  „  ,  an 

;  a  sy  ccmput  at  ions  show  that  e  A  -  ae  B  and  e  B  avf  :  r . 


e"  :  c  LA(-~,0; 


IR  )  where  g(  i  ) 


0.  The  main  result  of  this  gap  :  : 


Th'-trer;  ■  Suppose  that  (A)  is  satisfied.  Then 


1  ir  TN  ( r  )pN<  =  T(t  )<(>  ,  t  >  0,  <|>  e  Z  , 


it  :  <■  i  ng  uniform  on  bounded  intervals. 


The  resu  It  immediately  follows  from  t  Trot ter-Kato-The ore: 

a.  for  ir..  tar.se  contained  in  [1^1.  The  assumptions  of  this 
tr.eorem  will  be  verified  in  the  following  subsections. 


-*» k“«  .v. 7.S no,. .  v«  ^.  u  .  w. , 


*  .  V  *  •  •  v  V  *  V  MBA 


.1.  stability  of  the  scheme 

N  t- 

In  order  to  prove  |T  (t)|  <  Me“  ,  t  _>  0,  M  -  1  ,2  , 


show  that 


<  u!*||,  *  e  ZN,  N  =  1,2,..., 


where  u  is  independent  of  N  and  <•  •>.  rest.  I • 

g 

inner  product  reap.  correspond! ng  norm  on  Z  for  a  we. 

g  such  that  |*|“  and  |*L  are  equivalent,  |  <t>  I  ?  - 
C  g  Z  F 

j  ;T(i),c‘(T)i  di,  We  choose 


F(0  ;  F'O  +  Ip  j  1'or  T  6  (-hk,-hk-l  }  5 

J  -  K 

a  ;  l,...,p+l  (h  =  °°).  It  is  convenient  to  put  a,  = 

P  +  l  ,  k 

k  =  l,...,p.  Since  ( <t  °  ,  1  )  d> 1 )  e  ZN  for  <p  =  U°,<r 

O 

fjc  immediately  get 

1  NJ  p-_p*  i  T*  1 

<  6  x  >  $  5  2  =  <6  x,(l  +  s_k)<»  >  -  =  a.xp  (-h.-0), 
E  6 

\  1 

j  r  0 , .  .  .  ,p .  Using  this  formula  we  get  for  <p  e  Z  " 


<AN0,<t>^  =  [  A  «*.°  +  l  (A.  +  AnB.  )p  (-h  .  ) 
e  U  j  =  a  J  0  J  j 

+  /  (  A  (  T  )+  AqB(  t  )  )  <t>  ^  (  T  )dT 


+  J  g(T)(T7  ^ 1  )T<(|1  (  T  )di 


0  1  Pi  o 

+  a  U  -<>  (0)+  l  B .  <f>  (-h .  )+  / 

j=l  J  J  -CO 

+  I  Ua(-h,  )-*1(-h.-0)  ]T«1(-h.-0) 


r.ere  and  m  the  following  a^+i  =  1 .  We  shall  need  the  estimates 

0  , 
i  J  U(  T  )  +  A  B(  T  )  )4>  (  T  )dT  I 

—  Oo 

0 

<  ;  |  A  (  T  )  +  A  B  (  T  )  !  - —7p  g(T)1/2  U1(t)  'dT  (5.- 

U  g(T)1^ 

<  I  A  +  AnB  I  ~  <  ]  A  +  A  B  |  „  U1  !  _ 


0  '2  |V‘  r  2 
A-i.,  / 

1'  g  g 


0  T  2  ,y  1  2 
Ll/g  E 


J  B(T)«1(TjdTijc  |  B  j  U1!  2. 

■®  L.  ,  L- 

1/g  g 


Furt hermore 


/  g(T)(^7  )1  (t  )<!>1  (t  )dx  =  |  •  /  g(t)  U1(i)|2d- 


1  _ir  8  ( h+  t  )  d  i.l,  .  i2 ,  1  ?  _irj-l  d  (xl,  .,2, 

-  p  J  e  jTT  I  ^  ( x )  |  dT  +  -  l  a  jJ  — U  (  t  )  |  dx 

—  j=l  J-h.  Q 

J 

=  ilch-h-0)!2  -  I  ■/ 


( T  )  I  dT 


+  ^  a,  !  i^(-h  -0)  I2  -  1 4.1  (-h .  )  i  2  1  . 


"  j=l  J 

From  this  we  get 


J2  +  ^ 


|  .|a  -  ^(-hr0)|2 


+  t  j^VrV'^'V1 


+  |  1  «  UVh  -  -0 )  |  2  -  \  l  a  UVh  -0)| 

j  =1  J  J  j=lJ  J 

1|  Al/  H  n>|2  6  -1?  B(h  +  T  )  1  .  |  2 

+  -]4>  (-h-0)|  -  ^  J  e  1 4>  (t)|  dT 


-  32  - 


1  ,  1  i  p 

^  a.  i*  (0)  “ 
£  1 


f  "  ♦‘(-hj-0)' 


I  !  A,i.Vh.>ij  - 1  ‘/ 
j  =  l  J  J 


(  T  )  I  d  1 


For  I,  we  get  (also  using  (5-3)) 


1 

2  "l 


0  p  1  0  1  ? 
[  U  +  l  B  .  4>  ( -h  .  )  +  J  B(t)«  (x)dT  ' 

j=l  3  J 


_  U°  _  ha1  I2  _  |^(o)  |2 


❖  -  D  ❖  ‘ 


1  2a  (1  +  i  B  | 2  p  )|<t»|  l  +  a.  |  \  B.*1(-h.)|2 


Jl/g 


1  'j=l  J 


1  |  0 

-  —  !o 


n  ,  1  i  2  1.  l._v,2 

D<J  |  -  -^1  <p  (0)|  . 


Using  (5*2)  we  obtain  for  I  : 


X1  -  '  A0 


0 , 2 


*  U  *  A0B| 

1/g 


4°l  I*1! 


g 


1  '/V  \  |A.  *  A  B .  F  ♦  e  ! . 

J  0  J  j j 


4e 


I  u1(-h.)|2. 


By  hypothesis  (A)  we  can  choose  e  so  small  such  that 


J=1 

Alltogether  we  have 


B.*1(-h.  )  |2  -  A  \  A  .U1(-h  .)|2+e  I  U1(-h.  )  i2 
J  J  j  j  =  i  J 


.  N 

<A  $  > 


g 


(5 


-  (  I  A0i  +  !A"A0BiT2  +  ^liAi  +  A0Bi|2  +  2al(1+lBl?2  )}l* 


♦  e  zN. 


Jl/g 


.1=1 


Jl/g 


for  all 


The  proof  of  ( o .  1  )  in  case  of  ( A  .  18 )  is  quite  analogous.  Affair, 
one  uses  the  weighting  function  g  (instead  of  g(x)  =  e^r,  x  _< 
in  t  his  case). 


?.  Consistency  of  the  scheme 
I r.  this  section  we  Drove 


p  -*•  4'  as  M  -*•  °°  for  all  <t>  e  Z 


A  p  <p  -*■  A <t>  as  N  -*■  00  for  all  4>  e  V 


(5.5'. 


(5-6) 


where  f  c  don  A  is  such  that  (AI-A)D  is  dense  in  Z  for  A 
sufficiently  large. 

By  Lemma  2.3,  d)  the  sets  DR,  k  =  1,2,...,  satisfy  Dk  c  dorr, 
and- are  dense  in  Z.  since  (AI-A)Dk  =  (XI-A)Dk  and  =  dorr  A 

we  immediately  get 

(U-A)Dk  =  6k_r 

Therefore  the  sets  D  ,  k  -  1,2,...,  are  appropriate  candidates  f 

K 

v. 

Let  p ^ <p  =  (<P°,<pN)  for  4>  =  ( <}>0  .  Then  using  Proposition 

(A.  11),  (A.  12)  and  <t>  e  dom  A  we  get  for  4>  e  D. 


(A  p% 


Atf)°  =  [  (A.+A  B. )(«N(-h.  )  -  )) 

i -l  1  u  1  1 

0  N  1 

+  /  (  A(t  )  +  A  B(  t  )  )  (  <t  (  T  )  —  <t>  (l))dT, 


*  J  m  .  -  M  •  •  ^  f  4  VV,  -  V*  V  -  *_r  - "  *v  •-**. 


V  «:■  -  Ai ) 


i  u+  .  i; 


,N,  . i 


77  'f  _  $  +  5  U  4>  ( 0 )  -  $  (0) 

a  i  u 


1  B .  (  y  ( -h  .  )  -  y  ( -h  .  )  )-  /  B  '  T  )  (  4) 1  (  +  )  -  y 

i  =  l  1  1  1 


+  l  6'.  ( 4>  (~h.  )  -  y  (-h.  )  )  . 
i  =  1  1  1  1 


Therefore  by  (4.6)  and  (4.9) 


H+l 


fl"ph  -  -  (!fltA0Bi,2  *  -T7? 

Ll/g  rl 


[  \  I  1  N  j 

B  |  _  )  |  e  -  y  , 

T  ^  ]  c 

1/E  t 


+  "yuy(o)  -  y(0)i 


1/2 


i 

P-1  N+l  ,,n  ,  s  .  N 


+  .1  (  I  VA0Bil+£172(  •  Bi  I  +1 ) ) !  ^  (-h±  )-«>  (■ 


+  (!VVpl+1T72]Bpl  +  ^)1/2)U1(-h)-^( 


P-1 


+  (N  +  l)  l  U1(-hi)  -  *N(-h.-0)| 


1=1  ri+l 


+  (^)i/2iy(-h)-y(-h-o)i  +  i^(y-y)i  2 

L 


For  a  function  i>  e  L2  ( -<*  ,0  ;  IRn  )  we  introduce  ip .  -  4/|[-h.,-h_. 

.  e  i  n  1  i 

l  =  1,.  .  .  ,Pi  and  +  1  ~  Let  tt  i  =  resp. 

y  „  be- — -ii -  ^ 


.  wv.  the  orthogonal  projections  L  (-h.,-h.  ;  ]R  )  -+  Y.  resp 

-1  +  1  M  H  1  1  1  M  1m 


L^(-“»“h;  1R  )  -+  Y  ..  Furthermore,  we  denote  by  c>  resp.  o 
E  F  ^  n 

the  orthogonal  projections  L  (-1,1;  1R  )  -*•  span(P  I  , .  .  .  ,PMI )  resp 


LW(0,°°;1R  )  -*■  span(LQI , .  .  .  ,LNI  )  (recall  w(t)  =  e_t  ) . 


N‘ 


c  •  N  1 

Since  it  .  <j> . 

l  l 


.N 


i  -  l,...,p+l,  we  have  to  prove  that  as  N 


N  ,  1 


n  i  y  -  y  y  i 


TI  .  f  .  |  „ 

i  i  1 T  2 ,  ,  ,  mn  . 

L  (-hi  ,-hi_1  ;  TR  ) 


0,  i  1,...  ,p 


(5-7  ) 


n  i  y 

1  n  + 


p+  1 


p+i  p+i 1  L2(_a)j_h . 

s 


-  0 


(5-8) 


"i  +  l  i  -hi-0  '  ■*  0  >  1  = 


-  (i/V  )(-hi)  |  -  o,  i  = 


0  ,  i  =  1  , .  .  .  ,p  ,  (' 


d  .  kl  N  .  1  .  1 1 

n  '■  v  ]  77  i  ^  i  ^  1  p  n 

‘  1  Ld(-hi,-hi_1;  lRn) 


--(C1  -  7Tf'J  (J1  )| 

aT  P  +  1  P+1  P+1  Lg( >-h ;  lRn ) 

i  N  N 

for  4>  •  D.,  in  order  to  establish  |A  p  < 


It  is  easy  to  see  that 


?:  „  N  -1  .  .  N  .  n  '  N  -1 

i  -ic  °i  >  i  =  and  np  +  1  -  Gp  +  1o  0  p  +  1.  (d 


-*•  0  as  N 


Moreovo  r , 

-  / r'  i  \  1  /  2  „  |0  |  .  /  lxl/2 

i  ^ ;  i  (  2  '  5  1  ~  1  J  *  •  *  )P  3  ^p+1  ~  6  5 

( ! 

I  b  .  i  '  ( ^  )  ,1  1  ,  .  .  •  , p ,  I  + 1  1  "  8 

i  F 

For  •  l/  (-» ,0  ;  IRn )  we  put  x-  =  ©  •  1 4>  ^  ,  i  =  l,...,p  +  i. 

g  ill 

Lemma  5  •  2  .  a )  Let  $  e  D,  .  Then  x  •  c  (-1  ,1 ;  ]Rn )  ,  i  =  1  , 

_  r  k  _  1 

and  1  (  (<t1  )  ^  )  ,  j  z  0,...,k.  Moreover 

.1  1 

U|  -  =  (f)\((xi-oNxi)(j)),  j  =  o,...,k. 

i 


Let  <p  f  D  .  Then  x  ,  f  8, 

2k  p  +  1  k 

0  ,2  k.  Moreover, 


8v  and  =  (-  4  )jo'L  (u4 


8  p+1  "  rp  + 


0  , .  .  .  ,  2  k .  If  in  addition  <}>  6  E>2k+1  then  als°  xn  +  l  e  Bk‘ 


vp  +  l  k 


c )  Let  4>  e  Z  be  such  that  e  C ( - °° ,0  ;  IRn )  .  Then  e  C  ( —  1  ,1  ;  IF. 

i  =  l,...,p,  xp  +  1  e  C(0,°°;IRn)  and 

i  -  1 j  •  •  *  )P“^  j 

-  (^J)(-h.+0)  =  x±  (-1 )  -  (oNX;L)(-1),  i  =  l,...,p, 
4p*i('h)  -  (’p*-/p*i)(-h-0)  =  xPti (0)  -  <»NVi><°>- 

Prop f .  a)  is  an  easy  consequence  of  the  linearity  of  the  functions 

and  of  (5.13).  Similarly  we  get  the  formulas  for  the  derivatives 

of  -  and  of  the  error  <j>1  .  -  .<P'L  .  under  b).  c)  is  trivial. 

P  +  1  p+1  p+1  P  +  l^ 

It  remains  to  prove  x  .  e  B  for  <P  e  D01,  (and  x  „  e  8, 

p+l  k  kK  p+1  k 

for  4>  €  D2k  +  1).  Using  <P1  e  C2  k_1  ( -»  ,0  ;  IRn )  ,  (<t>a  )  (2k_1  ^  locally 
absolutely  continuous  on  (-“,0],  ( <J>p  +  1 )  ^  ^  e  L2  ( ,0  ;  IRn )  , 

j  =  0 , ...  ,2k,  we  get  xp  +  1  e  C2k-1 (0 IRn) ,  Xp2^_1 ^  locally 

absolutely  continuous  on  10,°°)  and  xp^  6  L2  ( 0 , 00 ;  IRn )  , 

j  =  0,...,2k.  For  m  =  0,1,2,...  and  j  =  0,...,2k  we  have 

(  t  +h  )m  ( <t>p  +  1  )  ^ ^  e  L2  (-°°,-h  ;  ]Rn)  by  Lemma  2.3,b).  Therefore 

tmx^j]  =  ('6)6Gp!l((T  +  h)m(4p  +  i)(j))  €  L2(0,~;]Rn).  By  Lemma  2.3,c) 
we ■ have 


lim  eBT/2(T  +  h)m(^  V)  =  0 

T  -4*  —  oo  P  ^ 


for  m  =  0,1,2,...  and  j  =  "0  , .  .  .  ,2k-l .  This  and  e~t/2tmx ^ ? (t ) 
<-6,”" L6h/:’o-t11(eB,/2(T*h)m(*p.1)(j))(t)  show  that  P 

lim  e  t/2tmxp‘)(t)  =  0.  Thi  ;  proves  x  .  e  B  .  The  proof  for 
t-co  P+1  p  +  1  k 

Xp  +  1  6  ^k  :''n  case  ^  e  D2k+1  ^ s  analogous  ■ 


2  and  Lemma 


Using  Theorem  3.1,  a)  for  k  =  2  and  Lemma  5.2,  a)  we  see  that 
( o .  7  )  id  satisfied  if  <*>  e  D_  ,  the  rate  of  convergence  being  r). 

<-  i  . 

Similarly  we  obtain  (5.8)  (by  Theorem  3-3  with  k  -  2  and 

1 

Lemma  5-2,  b))  for  $  e  D ,  with  rate  ^ ,  (5-9)  (by  Theorem  3.1,  c; 


N 


with  k  -  2  and  Lemma  5-2,  a)  and  c))  for  $  <  with  rate 

— ]  /  -  ,  (5-10)  (by  Theorem  3-5,  c)  with  k  =  1  and  Lemma  5-2,  b) 

and  c))  for  $  e  D,  with  rate  "  ,  (5.11)  (by  Theorem  3.1,  t) 

5  N  „  a 

with  k  =  2  and  Lemma  5*2,  a))  for  $  e  with  rate  — ana 


finally  (5.12)  (by  Theorem  3*5,  a)  with  a  - 

1 

Lemma  5.2,  b))  for  $  e  D„  with  rate  .  , 

■>  lj1/V 

shown  that 


,1/2 

1  ,  k  =  1  and 
).  Alltogether  we  have 


KT  r  4 

A‘'p  "4>  -  A*  =  0  ( -  1  0  )  for  e  D,, 

i.e.  (5.6)  is  established  with  V  -  . 

Condition  (5-5)  is  an  immediate  consequence  of  Lemma  5-2  a)  ar.d 
b)  and  the  completeness  of  the  Legendre  polynomials  in  L  (-1,1;  IF.) 

p 

resp.  the  Laguerre  polynomials  in  L^(0,®;]R). 

Therefore  all  assumptions  of  the  Trotter-Kato  theorem  in  Ilk) 
are  verified  and  the  proof  of  Theorem  5.1  is  finished. 


Remark .  If  B  =  0  then  <j>  e  can  be  replaced  by  <j>  e  D^, 

because  in  this  case  the  factor  N  is  not  present  in  (5.8) 

(and  in  ( 5  •  7  )  )  and  there  fore  <j>  e  is  sufficient  for  (5-8).  Of 
course,  the  smoothness  requirements  on  $  can  be  relaxed  if  one 
uses  interpolation  spaces  in  order  to  get  the  estimates  of  Section  3 
also  for  fractional  k. 

In  case  of  the  scheme  given  by  ( 4 . 1 8 )  we  get  for  <j>  e  D, 

K 


.  N  N  n  , 
A  p  $  -  A  ] 


i  ( :  A.A„B| 


(^V/2|Bi 


J!/P 


Jl/g 


-  0 


L‘ 


/  ,*  v  v  .•  .*  _•  . 

it  W*  WL  -V-*>.  V*  V.  d.  .» 


*  (!!l!)i'iV(o) -♦"«».  *  (fill)1'*  1  !b  I  |»1(-h, ) 

L  6  i=l  1  1 

-  *Vh.)i  +  l  I  A.  +  A0B.|  uVh.)  -  ♦N(-hi )  |  . 

J  ^ 

Proceeding  in  a  similar  way  as  above  we  get 

I  -  Al  |  =  0 (  -T-/V-)  for  <j>  e  D  . 

6  N  *  5 

The  reason  for  the  stronger  smoothness  requirement  in  this  case  i 

that  for  terms  like  ^  j  <d  (-bn  )  -  I***  ( -h .  )  |  ,  i  =  we  have 

to  use  part  b)  of  Theorem  3-5  instead  of  part  c).  If  B.  =  0, 

—  ~  J 

i  =  1  , .  .  .  ,p,  then  we  can  replace  D,.  by  D,. 

o  3 


5.3.  Approximation  of  the  nonhomogeneous  problem 

Since  (1.1)  is  linear  we  only  need  to  consider  the  case  <1=0 
and  f  i  0 . 


Proposition  3-3-  Let  z‘  (t)  be  the  solution  of  ( A . 1 3 )  and  let 
z(t)  =  (y(t),x  ),  x(t),  y(t)  being  the  solution  of  (1.1)  with 
1=0.  Then 

lim  2}'  (t  )  =  z(t  ) 

N-*-® 

uniformly  for  t  e  [0,tl  and  uniformly  for  f  in  bounded  sets  of 
L1  (0  ,t  ;  IRn)  ,  t  >  0  . 

Proof .  The  proof  is  analogous  to  the  proof  of  the  corresponding 
theorem  in  12],  using  the  variation  of  constants  formula  (2.6). 


5.^.  A  special  case 

The  scheme  presented  in  this  paper  has  the  remarkable  property 

to  give  the  exact  solution  in  specia]  cases.  Consider  (1.1)  with 

A.  =  B.  =  0,  j  =  l,...,p,  i.e.  we  have 
J  J 


D(xt)  =  x(t)  -  J  B  (  T  )  X  (  t  +  T  )  d  T  , 

0  -CO 

0  (5-!5 

L(x  )  =  A^x(t)  +  J  A(t )x(t+ t )dT . 

—  oo 

In  this  case,  the  schemes  defined  by  (A. 11),  (A. 12)  and  by  ( A . 1 8 ) 
coincide.  Furthermore,  the  t-method  as  described  in  [  7  ]  also 
yields  the  same  scheme.  We  put 


i 


i 


9 


a(x)  =  A(x)e  ^T,  b(x)  =  B(x)e  ,  x  _<  0 .  (5. It' 

Note,  that  (2.1)  is  equivalent  to  a,b  e  l?  ( -°°  ,0  ;  lFnx  n )  . 

£ 

Proposition  5 ♦ A •  In  addition  to  (5-15)  assume  that 
(i)  a,b  are  polynomials  of  degree  _<  m 

and 

1 

( I i )  is  a  polynomial  of  degree  _<  m. 

Let 

xN(t)  =  w!L(t)  +  l  B..w^(t),  t  >  0,  N  =  1,2,...,  (5.17) 

JU  j=0  J 

where  w\t)  =  col(w^(t  )  ,w^(t  )  ,.  .  .  ,wJJ(t  ))  is  the  solution  of  (A. 17) 
Then 

N 

x  (t)  =  x(t),  t  _>  0  ,  N  =  m,m+l,...  . 


Prop f .  Since  for 
the  columns  of  a 
for  j  _>  m+1 


j  _>  m+1  the  polynomials 
and  b  in 


L^(  — ,0;  IRn ) 


e  . 

J 

we  have 


in  (A. 19)  and  (A. 20).  Thus  for  any 

N 

(m+1 )n-dimensional  subspace  of  z  spanned  by  e 


are  orthogonal  tc 
=  0 

J  J 
N 


a  . 

J 

>  m 


=  3  • 
J 

the 


00 ,e0  ’  *  '  * 


m 


1  s 


invariant  with  respect  to  the  system  (A. 17).  Since  by  (ii)  we  have 
0  1  - 

(<J>  ,<(>  )  e  span{eQ0  ,eQ  , .  .  .  ,em  }  ,  the  coordinates 


0  ,  1 
of  the  solution 


N,.  . 
w  (t ) 


m 


N  .  N, 

woo  ^  ,w0 


M 


t  )  ,  •  *  •  ^  ^ 


of  (A. 17)  do  not  vary  with  N,  N  >  m. 


Let 


kN 


(t  ,  t  )  =  l  e  .  (  t  )w.  (t  )  ,  t  _>  0  , 
J'=0, 


to  Theorem  5-1  andJ "  Proposit ion  5.3 


<  0.  Then  according 


lim  wJJ  (t )  =  y(t  ) ,  lim  <J>N(t)  =  x.  in  L2 
N-*-“  uu  N->“  g 


(5.1“ 


uniformly  for  t  in  bounded  intervals.  Using  the  definition  of 
.  and  (i)  we  see  that 


N  N  ,  N 

x  (t)  -  Wq  (t)+  /  B(  t  )  ?  (t ,  t  )di 

—  CO 


t  ^  0)  N  :  m,m+l,...  .  This  shows  that 


xN(t)  =  xm(t),  t  _>  0,  N  _>  m •  (5.2:: 

From  (5.18),  (5-19)  and  (5.20)  we  obtain 

m  N  0 

x  (t)  =  lim  x  (t)  =  y(t)  +  /  B(t)x  (Odr  =  x(t) 

N-*’00  —00  ^ 

uniformly  for  t  in  bounded  intervals  ■ 


If  assumptions  (1)  and  (li)  of  Proposition  5.4  are  not 

M 

we  can  give  an  estimate  for  x(t)  -  x  (t). 


satisfied 


Proposition  5-5-  Consider  (l.l)  with  (5.15).  Let  *N  be  the 
orthogonal  projection  L|(-»,0;]Rn)  ■*  yN  =  span(eQ  , . . . ,eN)  and  put 


N 

a 


N 


11  a, 


Nl 

n  b. 


N  =  1,2,.... 


Then  for  any  t  >  0  there  exists  a  constant  c  not  dependent 
on  N  such  that 


I  x  ( t  )  -  XN(t)|  JC  c(  U1-  rrVl  +  |a-aN|  - 


b  -  b 


o 


) 


+ 


for  0  <  t  <  t,  N  =  1,2,...  ,  where  (t  )  i  .•  Fiver,  t  y  ( n  .  1':'  ■ 


Proof.  Let  T„,(t),  t  >  0  ,  be  the  solution  semigroup  generated 

-  N  — 

by  the  solutions  of  (1.1)  with  f  :  0,  A.  -  B.  -  0,  j  =  1  ,  .  .  .  ,p  , 

6tN  g  t  n  J  *1 

and  A,  B  replaced  by  e  a  ( t ) ,  e  b  (t),  respectively.  For  the 

same  equation  with  f  arbitrary  we  denote  the  solutions  correspo 

0  a  jj  0  N  1 

to  initial  data  <f>  =  ( $  ,  <t>  )  and  p'$  -  ( <t>  ,  *  )  by  x.(t),  yr 

and  x„(t),  y^(t),  respectively.  By  Proposition  2 . h 

t 

(y  (t),(x  )  )  =  T..(t)<*>  +  J  T  (t-s ) ( f (s ) ,0 )ds ,  t  >  0, 


(y  (t),(x  )  )  =  T  ( t  )  p  4>  +  J  T  (t-s ) ( f (s)  ,0)ds,  t  >  0. 


Proposition  5-^  implies 


x‘  ( t  )  =  x^ ( t  )  for  t  >  0, 


Using  the  second  equation  of  (1.1)  and  xM(t)  =  y.,(t)  + 
0  £t  N  -  N  h 

J  e  DTb  (  x  )x..(t  +  T  )dt  we  obtain 


x(t)  -  x  (t)  =  x(t)  -  x^(t) 


=  y(t)  -  yN(t)  +  J  e  Tb(T ) (x(t  +  i )  -  xN(t  +  i ) )dx 

—  OO 

0  Q 

+  j  eST(b(x)  -  b  ( T  )  )xN(t  +  T  )di 

—  ao 

0  fi  T 

=  y(t)  -  yN ( t )  +  J  e  b(x)(x(t+T)  -  xN(t+x))di 

—  oo 

0 

+  yN(t)  ■  yN(t)  +  /  e  u(t ) (xN(t+i )  -  xN(t+x ) )d- 


+  J  eBT(b(x)  -  bN(t )  )x  (t  +  T  )dT  ,  t  _>  0. 


!  x  ( t  )  -  x  (t)i  /2  max ( 1 , | B  |  ?  )t|(y(t),x  )  -  (y..!'t), 

IT.  . 

1/E  _  (5. 

+  |TN(t  )  U-pN$)  lz  }  +  |b-bN|  2  I  ( ^  N )  t  ■ ,  2 

L  '  L  ’ 

E  E 

for  t  >  0,  N  =  1,2,...  .  The  dissipat ivity  estimate  (5.^)  show 
that  there  exist  constants  M  _>  1  and  to  e  IR  such  that 

|TN(t) |  <  Meut,  t  >  0,  (5. 

for  all  N.  This  and  the  variation  of  constants  formula  imply 

|(x  )  |  <  MewtU|  +  M  /  eu,(t_s)  jf(s)|ds  (5- 

L  L  0 

g 

for  0  <  t  <  t  and  all  N.  From  Theorem  2.1,  c)  of  [91  we 
immediately  obtain 

l(y(t),xt)  -  (yN(t  )  ,(xN)t  |  z  <  c(|a-aN|  2  +  i  b  -  bN  :  0)  , 

L  L 

E  £ 

0  <  t  <  t,  where  c  is  not  dependent  on  N.  This  together  with 
(5-21)  -  (5.23)  implies  the  result  ■ 


Ir.  t  hie  section  we  discuss  some  numerical  examples  which 

demonstrate  the  feasability  of  our  scheme.  All  computation. 

were  performed  on  an  IBM  3081  at  Brown  University  usir.r  soft  war'. 

written  in  FORTRAN.  The  integration  of  the  system  (A.17)  of 

ordinary  differential  equations  was  carried  out  by  an  IMSL 

routine  (DVShK)  employing  the  Runge-Kutta-Verner  fifth  ana 

sixth  order  method.  The  coefficients  a.  and  6.  in  the  nr.'rix 

J  J 

H"  in  general  were  computed  using  lauss  quadrature  formulae  ;  q 


Ex  amp le  1 .  This  is  the  equation 

0 

x(t)  =  x(t)  -  j  (1-sini  )eTx(t  +  i  )di  ,  t  j>  0 , 

—  CO 

with  initial  conditions 

x(C )  =  1,  x(t)  -  0  for  t  <0. 

Because  of  the  special  initial  conditions  the  equation  is 
equivalent  to  the  Volterra  integro-dif ferential  equation 

t  _  , 

x  ( t  )  =  x  ( t )  -  J  (l  +  sin(t-t)e  T;x(i)di,  t  _>  0 
0 

x ( 0  )  =  1. 


(6.1  ) 


Differentiating  the  equation  in  (6.1)  we  see  that  the  solution  to 
(6.1)  also  satisfies  the  ordinary  differential  equation  (D  = 

(DU  +  2D5  +  2D2  +  D  +  1  )x ( t  )  =  0,  t  >  0, 

*  •  • ,  • •• 

x  ( 0  )  =  1,  x  ( 0  )  =  1,  x  (  0  )  =  0,  x  ( 0  )  =  -1. 

This  equation  was  used  in  order  to  compute  the  exact  solution  to 
problem  (6.1). 


Since  the  kernel  A(t)  -  (l-sint)e  is  oscillatory,  the  Ga: 

Laguerre  quadrature  formula  has  difficulties  to  yield  accurate 

values  for  the  a.'s  in  ( A . 1 9 ) .  However,  doing  the  same  compute' 

0 

as  in  the  proof  of  Theorem  3.5  >  c)  one  can  show  that 

OO  O 0 

I,  =  j  e  sin  wt  L,  (t)dt,  J .  =  J  e  cos  wt  L  (t)dt, 

K  0  K  K  0  K 

k  =  0,1,2... 


satisfy  the  recursion 


j  k  |  _  to 


-l)  fi  ) 
k-1  j 

“J  'Jk-J 


,  k  -  1,2,..., 


with  IQ  =  - =  - 2  *  Using  this  recursion  the  computatic 

1  +  OJ  1  +  u 

of -the  a.'s  posed  no  difficulties.  The  numerical  results  are  sh 

J  NO 

in  Table  6.1.  Note,  that  for  our  scheme  wqq(0)  =  4  for  all  N. 

Therefore  Table  6.1  does  not  contain  values  for  t  =  0. 


1.198396925 
1. 387642352 
1.5597599^ 

1.707896224 
1.826226289 
1.909918666 
1 .955135808 
1. 95  9053190 
l . 919884567 
2.0  |l. 836904312 


CPU 
( sec  ) 


1.198724502 
1. 389698823 
1.565189040 
1.717913392 
1 .841347388 
1 .929911297 
1.979088201 
1.985499876 
1 .946954731 
1 . 862474344 


0.029 


1.198671451 
1.389419452 
1 .564592523 
1.717066793 

1.840439101 

1.929179071 

1.978751749 

1.985714422 

1.947788385 

1.863905575 


0.054 


Table  6.1 


w32,t, 

"00vt' 

x(t  ) 

1.198669250 

1 .19866  9247 

1 . 389413284 

1 . 3894:3287 

1 .564588332 

1 .5645  8  :  341 

1.717070745 

1.71707:745 

1 . 840451863 

1 . 84045  :  65C 

1 . 929195874 

1.92919566: 

1.978765590 

1 . 97  87  655  69 

1.985719468 

1 .985716482 

1 . 947782026 

1 . 947762047 

1 . 863888825 

1.86386  6 83-9 

0.126 


•  this  ex  amp  le  the  assumptions  of  Proposition  5o  are  sat  i 

have  w  '_\  ( t  )  =  x^‘(t),  t  >  0,  because  b  =  0.  Observing  e 

uu  — 

obtain  from  Proposition  5-5  the  estimate 


x(t)  -  WQQ(t)|  £  const. |a  -  aN 


0  <  t  <  t  . 


It  is  easy  to  see  that  a(-x  )  e  8.  for  all  k  =  1,2,...  (note 

k  M  ^  v-  c-  f- 


„  .  ...  const 

t  =  1).  Therefore  according  to  Theorem  3-3  |  a  -  a  i  ~  L  - ~ 

L  ( M+ 1  )  *' 

all  k  -  1.2,...  (of  course  with  const,  depending  on  E  k)  which 
means  infinite  order  convergence  of  w^Ct)  -*■  x(t)  uniformly  o: 
compact  intervals.  This  is  reflected  by  Table  6.2  where  we  show 
Afi  =  max  |x(0.2i)  -  w!L(0.2i)|  for  N  =  4,8,16,32. 


g  k)  which 


i  =  0  , . . . ,10 


|x(0.2i)  -  wQQ(0.2i)|  for  N  =  4,8,16,32. 


N 

N 

A 

4 

0 .027897480 

8 

0.001414495 

16 

0.000016790 

32 

0.000000021 

Table  6.2 

The  next  two  examples  have  their  origin  in  the  dynamics  of 
structured  populations  (see  [13J).  Let  x  be  the  size  of  individua 
the  population,  0  <_  x  <_  S ,  S  being  the  maximal  size.  Then  a  simp! 
model  for  the  evolution  of  the  population  density  u(t,x)  is  give 
by 


u  (t,x)  +  (g(x)u(t  ,x)  )x  =  -  WqUO 

s 

u ( t  ,0  )  =  {  q(x)u(t  ,x)dx,  t  >  0, 
0 

u(0,x)  =  4>(x),  0  ^  x  <_  S . 


1 0  u  ( t  ,x),  t  _>  0  ,  0  _<  x  _<  S » 


Here  ^  >  0  is  a  mortality  rate  (assumed  to  be  constant),  g  is 
a  growth  rate  (assumed  to  be  positive  on  [0,S);  for  an  individual 


WWW 


-  46  - 


size  cha 
juried  to 
t:  i  a  1  size 


npes  accord inp 
be  nonne^at ive 
distribution 


to 


d  x 
dt 
and  es 


of  the 


-  g(x)),  q  is  a  fecur.i;  t.v 
sentially  bounded)  and 
popu lat ion . 


By  the  method  of  characteristics  one  can  show  that  the  Lirt 
rate  6 ( t )  =  u(t,0)  satisfies  the  Volterra  equation 

t 

B(t)  =  |  a(t-Os(C)dC  +  h(t),  t  >  0, 

0 

where 


aU) 

- 1  ~uo^ 

=  g(0)q(G  ( f )  )e  U  ,  £  >  0, 

hit) 

=  e 

U0^  ^  -1 

/  q(G  i(G(C)  +  t))«(Od 4 
0 

with 

GU) 

£ 

=  I 

0 

do 

- 

g(o  ) 

Assuming  that  h  is  locally  absolutely  continuous  or.  t  _> 
equation  (6.2)  is  of  type  (1.5). 


Example  2.  This  is  equation  (6.2)  with  =  0.15,  g(x)  -  b(S-x), 

b  =  0.0075,  S  =  60  and 


q(  t  ) 


-  2T_ 
? 

4S 


(“T^  +  St^ ) ,  0  <  T  <  60 . 


q  satisfies  q(0)  =  q(S)  =  q'(0)  =  q  '  (^  S  )  =  0,  q(^S  ) 
this  example 


27  -b£  2  ""  +  ^  ^  ’ 

a(C)  -  if  bS(l  -  e  L 


C  >  0, 


and 


yy •y  < 


lv> 


Vf. 


A  A' 


t  >  0 . 


+ 


-2  bt 


J  (1-|)  34>(  C)aC 
0  5 


Since  a(£)  is  a  linear  combination  of  e  J  ,  \ .  =  u  +  ib, 

J  0 

j  -  1,2,3,  8(t)  is  also  solution  of  an  ordinary  differential 

equation  (D  = 

dt 


(D+ A1  )(D+ >.p)(D+ A  )6(t)  =  ^  b3  S  6  ( t  )  , 
6(l)(0)  =  f(l)(0),  i  =  0,1,2. 


(£.- 


The  numerical  computations  were  carried  out  for  <Jj  =  l  or.  [C,S 
ar.d  the  results  are  listed  in  Table  6.3  (B(t)  was  obtained  by 
so lv ing  ( 6  .  4  ) ) 


\ 

4  ,  . 
“oo(t) 

“oo(t> 

' 

e(t ) 

:  i 

29.26615192 

29.26917027 

29.26917041 

29.269170^1 

i  ^ 

25. 38795761 

25.39133074 

25.39133081 

25.39133081 

1  3 

22.03843175 

22 .04086481 

22.04086408 

22.04086407 

1  4 

19.14858020 

19. 14960995 

19.14960983 

19.14960983 

i  5 

16.65708019 

16. 65672030 

16. 65672016 

16.65672016 

!  6 
t 

14 .50980845 

14.50831247 

14.50831234 

14.50831234 

7 

12.65928231 

12.65699712 

12.65699703 

12 . 65699703 

8 

9.68813725 

11.06133885 

11.06133881 

11.06133881 

9 

8.50035095 

9.68528995 

9.685289 96 

9.68528996 

10 

8.500350949 

8. 497625723 

8.497625775 

8. 497625775 

CPU 
( sec  ) 

0.021 

0.039 

0.087 

Table  6.3 

With  respect  to  the  rate  of  convergence  w^Ct)  -*■  x(t)  the  same 
remark  are  in  order  as  for  Example  1. 


- — — -  u 

:  0.0075,  0  -  60.  The  function  q 

Simple  comput  at.  ions  lead  to 


o 

is  the  same  as  for  Example  3. 


G  ( x  ) 


0  <  x  <  S, 


„-l 

U 


(y 


bSy 
1  +  by 


,  0 


y  <  ^ , 


and 


aU) 
h(  t  ) 


4 


27S 

4 


°  7.  ^  > 

(1+bC) 

|  (S-O  U  +  bt  (S-Q 
0  !S  +  bt(S-0]3 


i  0, 

0  (Ode , 


0. 


The  Lap lace-t rans form  of  a(£)  is  no  more  rational.  Therefore  the 
solution  6 (t )  of  (6.2)  in  this  case  does  not  satisfy  an  ordinary 
differential  equation.  The  results  of  our  numerical  computations 
are  shown  in  Table  6.4. 


s 

woo(t! 

"oh0 

1 6 , ,  s 

"oo(t) 

1 

1  3.  1  4  q  6  h  0  1  1 

I3.I556C69I 

13.15574424 

11 . 40999001 

11  .  4207 10  42 

11  .  42072946 

3 

9.95006974 

9.95575214 

9.95572216 

4 

8.72995239 

8.72962527 

8.72956630 

5 

7.71550920 

7.71OO8529 

7 .71002209 

6 

6. 87  4  8  0  2  8  5 

6.86594795 

6.86589901 

7 

6. 17887223 

6.16835637 

6.16833122 

8 

5 . 60209160 

5. 59145262 

5.59145294 

9 

5 . 12224059 

5.11264292 

5.11266498 

'10 

4.72038561 

4.71259466 

4.71263182 

CPU 
(sec  ) 

0.054 

O 

O 

0.129 

Table  6.4 


When  the  scheme  determined  by  (4.18)  is  applied  to  equations 
satisfying  (5.15)  the  only  quantities  to  be  stored  are  a.  and 


,  -  ‘"'ho  a;:  proximate  ly  linear  growth  of  CPU  time  observed 

U 

: 'at  ex  am:  irv  indicates  that  the  matrix  1A  1  is  not  stiff. 

The  next  two  examples  show  the  advantage  (as  far  as  rate  of 
convergence  is  concerned)  using  the  scheme  defined  by  (4.11), 
over  the  one  defined  by  (4.18)  when  (1.1)  involves  also  point 
delays  (see  also  the  remark  at  the  end  of  Section  5*2). 

Kxamn  1  e  1 .  m-lx ample  2.9  in  [10]).  This  is  the  retarded  problem 

0 

x(t  '  =  x(t-l)  +  2  J  seSx(t  +  s)ds ,  t  >  0, 

—  oo 

x  ( 0 }  -  .€  ,  x(t)  =  -t  ,  t  <  0. 


The  exact  solution  is  given  by 


x  ( 


1 


0  (t  )  =  --T  -  -t  -2  sin  t  + 


1  -2t 


x  v  t  )  - 


4  2' 

4 


4 


0  <  t  <  1  , 


4>(t)  +  y-  +  ^-(t-1)  -  ~r  cos(t-l)  -  ~  sin(t-l) 


+  (F(t-1  )cos  ( t  —  1  )  -  ^(t-1  )sin(t-l  )  -  -rgn  e  ^vt  1' 


.00 


20 


For  this  example  we  used  the  Legendre-Laguerr 


sc  he 


and  the 


Laguerre  scheme.  In  the  first  case  we  kept 


,.K 


.  -  span(e 


P 

->  r  5 


2,1 


only  increased  the  dimension  of  the  "Legendre"  subspace 
Y ^  =  span(e1  n,...,e1  M),  N  =  4,8,16,32.  The  results  are  liste 


1  ~r'“""'l,0**,**cl,N 

in  Table  6.5-  As  for  the  other  examples  A  is  the  maximum  of  the 


errors  at  the  meshpoints.  A  comparison  of  the  results  for  the 
Legendre-Laguerre  scheme  and  the  Laguerre-scheme  supports  the  rem 
at  the  end  of  Section  5-2. 
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Example  5  (Example  2.11  in  [10]  ).  This  is  the  neutral  •. 
equal i on 

d  0  s 

at  (x(t)  +  4  /  se  x(t  +  s)ds)  =  x(t-l),  t  _>  0 , 

—  CD 

x  ( 0 )  =  -4,  x  (t  )  =  1,  t  <  0  . 


The  true  solution  is  given  by 


x(t  ) 
x(t  ) 


*<t> 


v(t)  +  —  +  +  y^(t-l) 


0  <  t  < 
2 


27 


1  _-3(t-l)  4  N_  —  3(t  —  l ) 


+  9  6 


-  ^yCt-l  )e 


1  , 

t-1 

e 

1  <  t  <  2. 


The  computations  for  this  example  were  done  in  the  same  way  as 
for  the  previous  example.  With  respect  to  a  comparison  between 
the_  Legendre-Laguerre-scheme  and  the  Laguerre-scheme  the  same 
remarks  are  in  order.  The  CPU-times  and  the  errors  are  larger 
for  this  example  because  the  equation  is  of  neutral  type. 
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